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Abstract 

We discuss a neutrino mass model based on the S4 discrete symmetry where the 
symmetry breaking is triggered by the boundary conditions of the bulk right-handed 
neutrino in the fifth spacial dimension. The three generations of the left-handed 
lepton doublets and the right-handed neutrinos are assigned to be the triplets of 
54. The magnitudes of the lepton mixing angles, especially the reactor angle is 
related to the neutrino mass patterns and the model will be tested in future neutrino 
experiments, e.g., an early discovery of the reactor angle favors the normal hierarchy. 
For the inverted hierarchy, the lepton mixing is predicted to be almost the tri- 
bimaximal mixing. The size of the extra dimension has a connection to the possible 
mass spectrum; a small (large) volume corresponds to the normal (inverted) mass 
hierarchy. 



1 Introduction 



The origin of the fermion masses and mixing is one of the most compelling subjects in 
particle physics. In the past few decades, the neutrino oscillation experiments have re- 
vealed that the lepton mixing angles are quite different from the quark mixing matrix [1] . 
While the off-diagonal components of the quark mixing matrix are much less than unity, 
the lepton mixing angles are excellently modeled by the tri-bimaximal mixing |2J contain- 
ing two large and one vanishingly small angles. The observed generation structure has 
stimulated model building activities toward the theory beyond the Standard Model (SM). 

A phenomenological approach to explain the observed generation structure is to as- 
sume symmetry among three generations of fermion. Recently, discrete groups such as 
S3 [3], A4 [I], and S4 [5] have attracted much attention as feasible candidates for flavor 
symmetry in the lepton and quark sector. A common issue among these models is how 
to break flavor symmetry so that the observed data is naturally produced. Probably the 
most popular mechanism for flavor symmetry breaking is the vacuum expectation values 
(VEVs) of new scalar fields (flavons). Within this option, however, one is often forced 
to sacrifice the simplicity of the models due to multitude of flavons needed to construct 
desirable mass matrices. 

Another mechanism available in the literature is to utilize twisted boundary condi- 
tions [6] for bulk fields in extra spacial dimensions. This idea is explored within the see- 
saw mechanism in five dimensions where the S3 permutation symmetry is broken by the 
boundary conditions for the bulk right-handed neutrinos [7J. In this scheme, symmetry- 
breaking effects are calculable and the tri-bimaximal mixing is obtained as a consequence 
of symmetry breaking prompted by many possible boundary conditions. Furthermore, 
extra dimensions help VEVs of the flavons to be aligned so that viable neutrino mass 
matrices are obtained j8]. Since theories with extra dimensions have attracted great at- 
tention as a feasible paradigm to understand the gauge hierarchy problem [9j [TP] , the 
interplay between flavor symmetry and extra dimensions, it is an intriguing subject to 
investigate. 

In this paper, we discuss the neutrino model based on the five- dimensional framework 
of Ref . [7] . The novel point in this work is the adoption of the S4 symmetry as an alter- 
native to 5*3. Although the 64 group also belongs to the symmetric group as 5*3, it differs 
from S3 in many ways, so that the whole picture of the model and the physical conse- 
quences are significantly different from what is found in the previous study. For example, 
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the S4 group has various irreducible representations; two singlets, one doublet, and two 
triplets. By identifying the three generations of fermions as triplets, the number of the 
Yukawa couplings is restricted to be small, so that the model prediction becomes rich and 
testable at future experiments. In addition, the variety of the irreducible representation 
makes it possible to construct realistic charged-lepton mass matrices. Another point is 
the order of S4 (4! = 24) which is much larger than S3 (3! = 6). With the 24 group 
elements in S4, the number of the theoretically possible boundary conditions for the bulk 
fermions is large (100 patterns) and the theory contains new types of the lepton mass 
matrices which are unavailable in the S3 model. 

This paper is organized as follows. In Section 2, we discuss the basic setup of the 
model; general boundary conditions for the bulk neutrinos, locations of the SM fields, the 
Kaluza-Klein (KK) expansion etc. In Section 3, we present two concrete models of the 
S4 breaking and discuss the predictions which can be used as markers in future neutrino 
experiments. Section 4 is devoted to summarizing the results. Appendices |A] and IB1 
summarize our convention for Lorentz spinors and some basic properties of the S4 group, 
respectively. 

2 Basic framework 

We consider the fifth spacial dimension y compactified on the S 1 /Z 2 orbifold with a radius 
R. There exist two fixed points at y = and nR. The four-component bulk fermions 
*&i(x, y) (i = 1,2,3) are introduced as three generations of the right-handed neutrinos and 
their chiral partners [H] . There are two boundary conditions linked to the two operations 
on the S 1 /Z 2 space; the translation T : y — )■ y + 2nR and the parity Z : y — >■ —y. The 
boundary conditions for the bulk fermions *&i(x, y) are written as 

%{x,-y) = Z ij ®7 l & j {x,y), $;(x, y + 2ttR) = T i3 ; ^(x, y), (2.1) 

where Zij and T^- are the matrices acting on the generation space. The parity and trans- 
lation imply that Z 2 = 1 and ZT = T~ l Z must be satisfied. Instead of the translation 
Tij in (12.11) . one may use another parity Z' = TZ to write the boundary conditions 

M x >-V) = z v^'y5%(x,y), V i (x,irR-y) = Z< j ®'y B V j (x,irR + y). (2.2) 

The parity Z-j is the reflection with respect to y = ttR. In this case, the two matrices 
must obey Z 2 = 1 and Z' 2 = 1 as consistency relations. 
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The SM fields including the left-handed lepton doublets = PlLi = ( ) and 
the Higgs field H are assumed to be localized at y = nR. This SM-field profile gives an 
example and the analysis below is applied to other cases in a similar way. The Lagrangian 
including the neutrino fields is given by 



^jV'Wj - ^(M^. + h.c.) 



1 / 



y/A V 



VitY^LjH + 9?K t )LjH + h.c. 6(y - «R) 



(2.3) 



where A stands for the fundamental scale of the theory. In this paper, we do not consider 
Lorentz- violating mass terms such as \I/r 5 \I/, \l/ c r 5 \l/ [12] and the bulk Dirac (kink) mass 
while they are often discussed in literature. After the electroweak symmetry breaking, the 
boundary interactions develop the neutrino Dirac masses; = Y n .v and m\j = Y£ v 
where v is the VEV of the Higgs field (H) = ( ° ) . The charge-conjugate spinor \l/ c is 
defined by \l/ c = T 3 r 1 \l/ . Our convention for the gamma matrices and Lorentz spinors 
are given in Appendix |A] 

The four-dimensional effective theory is described by the KK modes of the right- 
handed neutrinos. Under a set of boundary conditions, the bulk fermions ^(x, y) are 
expanded as 



/ \ 

1 Y,^kM^U x ) 



n=0 

oo 



\ n=0 



(2.4) 
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with the orthogonal systems Xrl(v)- ^ * s convenient to choose them to satisfy the nor- 



malization conditions J^ R dy 

mode is canonically normalized. By substituting the expansion into the five- dimensional 
Lagrangian and integrating it over the extra space, we have the four-dimensional effective 
Lagrangian 



m ] n 
AR,L X.R, L 



fimnfiij so that the kinetic term of each KK 



£ 4 = iN^a^d^N - - (N T e g> M N N + h.c.) 



(2.5) 
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where e is the antisymmetric tensor and 

tR 



M Kmn = I* dy xllH-dyM, 
Jo 

wR pttR 

WlTji/r„,n /\/f _ / J„, ,,raT; 



M Rmn = / dyx^Mxl, M Lmn = / dyxTMxl, 
Jo Jo 

M n = -±= X rH*R)™, M c n = -±= X n J{*R)™ c . (2.7) 



The zero modes are suitably subtracted according to the boundary conditions. The gen- 
eration indices in (12. 5p . (12. 6p . (12. 7p are suppressed for simplicity. We leave the entries 
connecting, different KK levels blank since these entries will be vanishing in the following 
discussion^. 

The mass spectrum of Majorana neutrinos is obtained by diagonalizing Mn- For 
Mp <C Mn, the seesaw mechanism is available and the Majorana mass matrix of the light 
neutrinos is approximately given by 

M v = —MpMjj 1 M D . (2.8) 

In what follows, we assume that the seesaw mechanism works. That is, the inverse of 
the compactification radius and/or the bulk Majorana scale are much larger than the 
boundary Dirac masses M n and M^. 



3 5*4 symmetry breaking and its consequences 

In this section, we introduce the S4 discrete symmetry and calculate its breaking effects 
on the neutrino mass matrix. We present S4 charge assignment and the possible boundary 
conditions allowed by the consistency conditions. Two particular boundary conditions are 
discussed in detail as viable examples. Finally we comment on possible constructions of 
the charged-lepton sector. 

3.1 Charge assignment and possible boundary conditions 

The irreducible representations of S4 are two singlets 1 and 1', one doublet 2, and two 
triplets 3 and 3' (see Appendix [B]). Suppose that the three generations of the bulk fermions 

*This does not hold in more general/different setups, e.g., on curved background or the boundary 
Majorana mass [15] . 
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Table 1: The possible boundary conditions in terms of the two parities Z and Z ' . The 
symbols a\, 0,%, a 3 , ■ • ■ , e 4 stand for the group elements (see Appendix [B]). The calligraphic 
characters represent the boundary conditions with which the theory becomes viable for 
neutrino phenomenology. Out of the 100 general possibilities, 48 patterns are useful. 



^i(x,y) and the lepton doublets Li(x) behave as triplet 3. The symmetry-invariant mass 
parameters in the five-dimensional Lagrangian ( 12. 31) are then written as 



Mi A = MSi 



1.1 ■ 



rriij = mSij, 



m c b, 



13 ■ 



(3.1) 



If symmetry is preserved, the neutrino mass matrix M v is also proportional to the iden- 
tity matrix, which is not consistent with the observations. In fact, the trivial boundary 
conditions Z = 1 and T = 1 lead to the neutrino mass matrix [13] 

•2 



1 \M\R rrr 
" y AR tanh(7r|M|i?) M* iy 



(3.2) 



Symmetry breaking is the key to obtain the flavor mixture and the mass splittings between 
three generations. 

It is convenient to specify the boundary conditions by Z and Z 1 with which the consis- 
tency relations are written as the parity conditions Z 2 = 1 and Z' 2 = 1. Tabled] shows the 
possible combinations of Z and Z'. In the S4 group, there are 10 elements satisfying the 
parity condition, so that there are total 100 possibilities in the table. The boundary condi- 
tions are classified into 14 categories according to their physical consequences. Out of the 
14 categories of the boundary condition, the three categories tagged by the calligraphic 
characters A, B, C prompt S4 breaking which is viable for neutrino phenomenology. The 
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subscripts are the serial numbers in each category. The conditions in the same category 
produce identical neutrino mass matrices up to the rotations by the group elements. For 
instance, the seesaw induced mass M v for the condition B\ is obtained by exchanging the 
second and third generation labels of £>2, and vice versa. Since such rotations are absorbed 
by appropriate S4 transformations of the left-handed neutrinos, physical implications in 
each category are equivalent. 

Within the three useful categories of the boundary conditions, the S4 symmetry is 
completely broken and the mass matrix at the low-energy acquires structure sufficiently 
rich to account for the observed neutrino masses and mixing. Since the neutrino mass 
matrix with the tri-bimaximal mixing is invariant under the S2 transformations, one may 
naively expect that 5*4 — > S2 breaking patterns are feasible for neutrino phenomenology. 
However, such categories produce too simple mass matrices to be realistic. For example, 
many S4 — > S2 breaking conditions predict that at least two neutrino masses are degen- 
erate. In what follows, we discuss physical implications of the three categories A, B and 
C in detail. 

3.2 Model I: The category A 

Let us first discuss the category A. Out of the 24 patterns in A, we focus on Ag. 



fi 




Z' = dt 




(3.3) 



as a concrete example which is convenient for presentation. The KK expansion (12. 4p 
which satisfies the boundary conditions is given by 
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(n > 1), (3.4) 



where V is the unitary matrix 



V 




(3.5) 



with to = e i2w / 3 . Three "zero modes" are absent in this boundary condition (here 
^2 ' ^°l 2 ' ^L a are taken as such absent fields). In order to satisfy the boundary condi- 
tion, nontrivial generation structure is necessary in the KK wave functions. 

By substituting these KK expansions into (12.71) and performing the seesaw opera- 
tion (I2.8p . one finds the Majorana mass matrix 
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(3.6) 



where c = cosh(27r|M|i?) and s = sinh(27r|M|i?). It is noticed that the 5*4 symmetry 
is completely broken in (13. 6p . However, in the limit |M|i2 1, the last term becomes 
negligible and M v recovers the Si symmetry (2o3 exchange). This is because the theory 
has two breaking sources at y = and y = ttR, and £4 is entirely broken only globally. In 
the large-size limit of extra dimension |M|i? 3> 1, the boundary condition at the distant 
brane (fi at y = 0) becomes irrelevant to physics at another boundary, whereas the local 
twisting (d\ at y = ttR) which respects the 2-3 exchange symmetry remains relevant. 

If the last term of (13.61) is vanishing, the Majorana mass matrix (13.61) is diagonalized 
by the tri-bimaximal mixing [2] 



Vtxi 
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(3.7) 



Besides the large-volume limit mentioned above, that is also realized by taking m c = 
which predicts the inverted hierarchy with m 3 = 0. For the normal hierarchy, both m^O 
and m c 7^ are necessary, so that the neutrino mixing is deviated from the tri-bimaximal 
form by presence of the last term. The deviation induces a nonzero reactor angle #13 
sufficiently large to be measurable in forthcoming experiments. 
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Let us examine the mass matrix (13.61) and its predictions in detail. It is useful to 
rewrite the matrix (13.61) as 



M„ 
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The neutrino masses are given by 
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If the charged-lepton mass matrix is diagonal, the lepton mixing 



matrix U is identified as the unitary matrix which diagonalizes (13. 8ft ; 
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where p = arg(m) + arg(M)/2, a = arg(m c ) — arg(M)/2 and 

tan26» = , ^— ox . 

s(l + r 2 ) 

The relevant mixing matrix elements are written as 
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The solar angle is robustly predicted to be sin #12 = 1/ v^3 up to the small corrections of 
0(9f 3 ). The mixing matrix U contains the phase parameters which cannot be removed by 
field redefinitions. That is, the boundary condition induces not only the S4 breaking but 
also CP violation. The mass formulas ( 13 .9p involve three effective parameters; |m| 2 /A, 
|M|i? and r. These parameters are determined if the three neutrino masses are regarded 
as fixed observables. 

The neutrino mass matrix ( 13. 8 h accommodates all possible mass patterns of neutrinos; 
normal, inverted and degenerate. The mass pattern is determined by the mass ratio r. 
In the region where r > 1 (r < 1), the normal (inverted) mass ordering is realized. In the 
case of r = 1, neutrino masses are degenerate since \m\\ = \ 777.3 1- 
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Figure 1: The reactor angle sin #13 as a function of mi + m2 + for the normal mass 
ordering in the Model I. 



At first, let us discuss the normal mass ordering (r > 1). From the mass formulas (13. 9p . 
it is seen that a hierarchical spectrum is realized in the regime that r > 1, \M\R <^ 1 
and \M\Rr > a/3/2tt. In such a regime, the ratio of the two mass squared differences is 
approximated as 

8x 3 V3 
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9 - 3x 2 - 8x 4 + 4xV3 + 4a; 2 ' 

where x = ir\M\Rr. By substituting the typical observed value lAm^l/Am^ - 
the left-hand side, one obtains x ~ 1.95. The reactor angle is then predicted as 
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(3.13) 
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(3.14) 



Interestingly, the reactor angle is predicted to be just below the current upper bound. 
The discovery of #13 is imminent if the normal hierarchy is realized in this model. It is 
seen in ( 13. 14|) that #13 decreases as the parameter r decreases, which means that #13 takes 
maximum value at the hierarchical limit and it is lowered as the spectrum is shifted to 
the degenerate pattern. 

Figure [1] shows the numerical plot of sin # 13 as a function of the total neutrino mass 
£ = X]i m «- A s seen i n (13 -9ft and (13. lip , the three parameters |m| 2 /A, |M|i2 and r are 
determined by regarding jAm^l, Am^ and S as inputs. That is, sin 6*13 is plotted as a 
function of £ with fixed values of IAtttJJ and Am^. In the numerical calculation, the 3cr 
ranges of the mass squared differences in Ref. [H] are used. The cutoff behavior of ( 13 . 141) 
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Figure 2: Correlation between sin 2 2023 and sin #13 for the normal mass hierarchy in the 



is seen at the lower end of E, i.e., at the hierarchical limit. In the hierarchical regime 
that y|Am|7f ^ S < 0.06 eV, 0.07 < sin #13 < 0.17, which is expected to be observed 
at T2K, Double Chooz, RENO and Daya Bay [Tjj]. As £ increases, that is, the three 
neutrino masses approach the degenerate pattern, sin #13 becomes smaller. This behavior 
is also understood in view of (13. 8p . It is seen that the three diagonal components get 
closer to each other in the regime |M|i? 3> 1 and r ~ 1. The off-diagonal elements are 
then suppressed by a large cosh(27r|M|i?) factor compared to the diagonal ones. 

The atmospheric angle #23 is, as seen in (I3.12p . correlated to #13. If CP is conserved, 
623 is deviated from the maximum in such a way that sin 023 ^ 1^/2 + sin 0i 3 /2. With 
a finite combination of the phase parameters p — a, however, the effect of #13 can be 
canceled and the maximal 023 is possible even with a nonzero value of #13. Figure [2] 
presents the correlation between sin 2 2023 and sin #13 where the shaded region shows the 
possible parameter space. The lower bound of #i 3 is predicted for each fixed value of 
sin 2 2#23 . In addition to the correlation between sin #13 and the mass pattern, this will 
be also a crucial test of this model if sin 2 2# 23 will be precisely determined at the T2K 
experiment. 

Next let us discuss the inverted mass ordering (r < 1). As seen in (13. 8ft and (13. 9p . m\ 
and ?7i2 are close to each other if |M|i? ^> 1. In this regime, tanh(7r|M|i?) ~ 1 — 1/s and 
the ratio of the neutrino mass squared differences is given by 
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By using this relation, the reactor angle sin #13 is written by r and the mass differences; 
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Even at the extremum r = l/VS, sin #13 = 0(10 3 ) which is far below the expected 
sensitivity in future experiments. After all, the reactor angle follows 



sin #13 < 0.003 



(3.17) 



with typical values of the mass differences. The lepton mixing matrix is thus almost the 
tri-bimaximal mixing in the case of the inverted mass hierarchy. 



3.3 Model H: The categories B and C 

Another interesting example is the categories B and C. These two categories lead to quite 
similar physical consequences. The predictions differ only in the solar angle 612 ■ Hence we 
regard B and C as a single model and discuss B in detail. We demonstrate the viability 
of the model with a concrete example B$: 
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where O is the orthogonal matrix 
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(3.23) 



Two "zero modes" are absent under this boundary condition (here ip^ , i\P L are taken as 
such absent fields). With the KK functions fl3.19p - fl3.22p . the neutrino mass matrix at 
low energy turns out 
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In this model, the solar angle 812 originates in the neutrino sector while the reactor and 
atmospheric angles are generated in the charged-lepton sector. A suitable charged-lepton 
mass matrix will be discussed in the next subsection. For the category C, the mass 
matrix (I3.24p is modified by the 45° rotation for the 1-2 sector. 

By diagonalizing the upper-left 2x2 submatrix in (I3.24p . one finds the mass eigenvalues 
and the mixing angle 
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This model accommodates only the inverted mass ordering since m 3 cannot be largest 
among the three eigenvalues. The two eigenvalues m± are identified as m<i = m 



+(-) 



and 



mi = ^i-(+) for r > 1 (r < 1) and degenerate for r = 1. The observed mass differences 
are not reproduced in the regime that |M|i? > 1/27T. In such a regime, the ratio of the 
mass differences is given by lAm^l/Am^ ~ — 1 for r > 1 and 2r 4 /s for r < 1, which 
contradict the data. 

A realistic mass pattern is obtained if the bulk Majorana scale is much smaller than 
the inverse of the compactification radius: \M\R <^ 1. In such a regime, the parameters 
r and \M\R are related as 
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Figure 3: The solar angle sin 2 9 12 as a function of m\ + mi + m% in the Model I. The 
horizontal lines shows the 3a allowed range of the solar angle [T4"] . 



where the plus (minus) sign in the right-hand side is for r > 1 (r < 1). By using (I3.26p . 
the magnitude of 6*12 is estimated as 
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tan2# 12 ~ 4 
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r 2 + 1 

The angle #12 is zero for r = 1 and it rapidly increases as r departs from unity due to the 
large factor lAmfjJ/Am^ ~ 30. That is, the solar angle is small at the degenerate limit 
while it becomes large as the spectrum approaches the inverted hierarchy. 

Figure |3] shows sin 2 9 12 as a function of the total neutrino mass E = J2i m i- The blue 
(shaded) region is for the category B and the pink (light-shaded) region is for C. The 
horizontal band is the 3a range of sin 2 # 12 [2j. The thickness of each plot corresponds to 
the possible 3a values of \Am\ 1 \//\m\ 1 in the same reference. 

In the case B, the solar angle is increased as the spectrum is shifted from the degenerate 
pattern to the hierarchical one. The 3a range of sin 2 9 12 = 0.27 — 0.38 is reproduced for 
O.lOeV < E < 0.12 eV, which range will be tested in the cosmological observation such 
as cosmic microwave background data from the Planck mission. While in the case C, the 
solar angle is increased as the total mass is increased. This opposite behavior is because 
of the 45° difference of the solar angle 9u- The predicted range of E is slightly wider than 
the case of B and shifted to a higher region: 0.11 eV < E < 0.15 eV. 

In this model, 9 13 and #23 comes from the charged-lepton sector and their magnitudes 
depend on the structure of the charged-lepton mass matrix. The charged-lepton mass 
matrices suitable for the two options of the neutrino sector (Model I and E ) are discussed 
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in the next subsection. 



3.4 The charged-lepton sector 

The charged-lepton sector resides on the SM boundary y = 7rR. The right-handed electron 
e R and the pair of muon and tauon (/j, R , t r ) are assigned to be 1 and 2, respectively. We 
also introduce new gauge singlet scalars 0j, which are assigned to be 3. These assignments 
are summarized in the table below: 





£r (vr,tr) (L e , Lp, L t ) 


H 


(01,02,03) 


S4 


12 3 


1 


3 



For the charged leptons, the 64 invariant Lagrangian is 

C = (L e 0i + L M 2 + L r 3 ) H* 

+ ~£ [ ('PR + ^r) L e0i + (^ 2 7Ar + ujt^) ^02 + (wpR + w 2 r^) L r 3 ] H* , (3.28) 

where Y s , Yd are the dimensionless Yukawa couplings. After the electroweak symmetry 
breaking, the charged-lepton mass matrix is given by 

(«i a 2 a 3 \ / \ 

+ vY d I ai Lu 2 a 2 wa 3 , (3.29) 
J \a± ua 2 io 2 a 3 J 

where = (0j)/A. The Hermitian matrix M\Mi becomes 

/(|y s | 2 + 2|r d |> 2 (|r s | 2 -|r,| 2 )« ia2 (|y s | 2 - \Y d \ 2 ) ai a 3 \ 

M}M t = v 2 (|y s | 2 -|r d | 2 ) ai a 2 (|F s | 2 + 2|F d | 2 )a 2 (|F s | 2 -|F d | 2 )a 2 « 3 . (3.30) 
V(in| 2 - |Ki| 2 ) « lCK3 (|n| 2 - > 2 a 3 (|r,| 2 + 2|F ci | 2 )a 2 / 

By assuming that a%v ~ m e , a 2 v ~ m M , a 3 u ~ m T and the Yukawa couplings are of order 
unity, one obtains the realistic charged-lepton masses and small mixing in (I3.30p . This 
case is suitable for Model I for the neutrino sector because the near tri-bimaximal mixing 
comes from the neutrino sector. Corrections from the charged-lepton mass matrix are at 
most 0{m e /m l j), 0(m e /m T ) and (9(m^/m r ) for the lepton mixing angles 812, #13 and 823, 
respectively. 

On the other hand, the neutrino sector of Model I provides only Qyi and needs (at 
least) a large 2-3 mixing arising from the charged-lepton mass matrix. This is possible if 
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ai\ — and a 2 = «3 = ct are realized; 

/O \ 

MtM, = ^ O (|F S | 2 + 2|F,| 2 ) a 2 (|F S | 2 - \Y d \ 2 ) a 2 . (3.31) 
\0 (|y s | 2 -|y,|> 2 (|F s | 2 + 2|F d | 2 )aV 

The charged-lepton masses and the mixing angle are 

7T 



m e 



= 0, mJ = 3|F d | 2 aV, m 2 = (2|F S | 2 + \Y d \ 2 ) a 2 v 2 , e 23 = J . (3.32) 



By taking Y d /Y s ~ sj2/?>m j j i /m T , one obtains the observed mass ratio m M /m T and the 
large 2-3 angle which accounts for the atmospheric neutrino oscillation. For the Yukawa 
coupling Y s of order unity, a = O(10~ 2 ). The nonzero electron mass is produced by 
holding a finite but small value of a±. 

Another interesting way to obtain the realistic electron mass is to take account of the 
higher-order corrections to the charged-lepton mass matrix while holding the vacuum- 
alignment aii = and a 2 = a 3 = a strictly. The Lagrangian with the dimension six 
operators is given by 

Y S A _ 

A£ = ~^e R ( L e 2 03 + £ M 0103 + £r0102) H* 

Y A _ _ _ _ _ 

+ [ (HR + T R ) L e 2 03 + (w 2 /!^ + LOTr) L^010 3 + (u/jLr + UJ 2 Tr) L r 0l0 2 ] H* . 

(3.33) 

After (pi developing the VEVs of {(pi) = (0,aA,aA), the charged-lepton mass matrix up 
to this order is written as 

/Y A a 2 Y s a Y s a \ 
M e = M E + AM, = v Y d A a 2 Y d u 2 a Y d ua . (3.34) 

\Y A a 2 Y d ua Y d u 2 a) 

t 



The Hermitian matrix Ml Me follows 



'(|F s A | 2 + 2|F f A | 2 )a 2 (Y A *Y s -Y A *Y d )a {Y A *Y S -Y A *Y d ) 
M\M t = a 2 v 2 [ (Y A Y* - Y A Y*) a \Y S \ 2 + 2\Y d \ 2 \Y S \ 2 - \Y d 

(Y A Y s *-Y A Y*)a \Y s \ 2 -\Y d \ 2 \Y S \ 2 + 2\Y C 



2 



(3.35) 



The charged-lepton masses and mixing are given by 

m\ ~ 2|r d A | 2 aV, m\ ~ 3\Y d \ 2 a 2 v 2 , m 2 ~ (2|r s | 2 + \Y d \ 2 ) a 2 v 2 , (3.36) 



y/2\Y A Y s * - 


Y A 


d | 


a 


2in: 


2 + 


Y d \ 


2 



7Y 



^12=0, ^ 'o^, S 2 I , ^23 = p (3-37) 
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where the 1-3 angle 6[ 3 is assumed to be small not to contradict the reactor bound. As 
in the case of the leading-order estimation, Y d /Y s ~ a/2/3 m M /m T and a = 0(1CT 2 ) for 
Y s ~ 1. The electron mass is in the right range if Y d A Y^. The 1-3 angle becomes 
0f 3 ~ l/V2(K s A /y s ) which will be observable if Y S A /Y S = C^ICT 1 ). 

A possible way to achieve the aligned vacuum = (0, aA, aA) is to assume that 
the scalar fields <$>{ propagate in the five-dimensional bulk and follow nontrivial boundary 
conditions [8]. By assigning the pseudotriplet 3' to 4>i or assuming the scalar potential 
to have additional Z 2 symmetry, it is possible to impose the twisted boundary condition 
<f)(y + 27r i?) = — G?20(y). While the first component (f>i cannot develop a constant VEV 
because of the antiperiodic boundary condition, the other two components 02 and 03 
acquire equal VEVs since the boundary condition respects S2 symmetry which is the 
exchange of the 2-3 component of the triplet fy. 

The bulk scalars <pi with nontrivial boundary conditions are also useful for realizing a 
hierarchal charged-lepton mass matrix suitable for the Model I. Suppose that the right- 
handed charged-leptons (e#, t r ) are 3 of S4 and also propagate in the bulk. 

The gauge singlet scalars (pi are also coupled to the neutrinos via higher-dimensional 
operators, e.g., ^>LHcj). For the Model I, the most stringent contribution comes from (<j) 3 ) 
component related to the tau mass: m T ~ -\/l^s| 2 + 2|Fd| 2 t>a3. However, a 3 = 0(1CT 2 ) if 
all the dimensionless couplings are of order unity. Thus the corrections from such higher- 
dimensional operators remain a few percent and the predictions discussed in Section 13.21 
and 13.31 are valid unless the Yukawa couplings are anomalously small. A similar discussion 
is also hold for the Model E. 

4 Summary 

In this paper, we have examined the S4 discrete group as a feasible candidate for the fam- 
ily symmetry which accounts for the leptonic generation structure. Phenomenologically 
viable S4 breaking is triggered by the boundary conditions of the right-handed neutrinos 
which reside on the five- dimensional spacetime. With the SM fields being trapped in 
the four-dimensional subspace, many boundary conditions produce realistic lepton mix- 
ing slightly deviated from the tri-bimaximal mixing in a controllable way. There are two 
types of plausible structures (Model I and E in Section 13.21 and 13.31 respectively) in the 
neutrino sector. 

In the Model I, the whole mixing matrix originates in the neutrino sector and the 
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charged- lepton mass matrix is diagonal in a certain S4 basis (see Appendix [Bj. While the 
solar component U e2 is robustly predicted to be 1/ a/3, the reactor angle #13 is correlated 
to the mass patterns. With the normal mass ordering, # 13 increases as the total neutrino 
mass decreases in such a way that sin^ > 0.07 if mi + m<i + 777,3 ^ 0.06 eV, which range 
will be probed at future experiments. Furthermore, the atmospheric angle #23 is deviated 
from the maximal value 45° in such a way that the deviation is ruled by the magnitude 
of #13. While the discovery of #13 is imminent for the normal mass hierarchy, the reactor 
angle is generally small (6 I 13 < 10 -3 ) for the inverted mass ordering; the mixing matrix is 
predicted to be almost tri-bimaximal form in the case of the inverted ordering. 

While in the Model E, the large solar angle (9 12 stems from the neutrino sector and the 
maximal atmospheric angle #23 emerges from the charged-lepton mass matrix. Unlike the 
Model I, the Model 1 accommodates only the inverted mass ordering. The solar angle Qyi is 
related to the total neutrino mass. Within 3a range of the current oscillation parameters, 
the total mass is predicted to be in a small window 0.10 eV < mi + m 2 + m 3 < 0.15 eV 
which can be probed by the Planck satellite. By the combined data of future oscillation 
experiments and cosmological observations, the models discussed here are not only tested 
but also distinguished each other in their predictions for the leptonic generation structure. 

Besides the predictions in the masses and the mixing angles discussed in this paper, 
there are other phenomenological issues such as collider signals [IE], leptogenesis [17] , 
low-energy CP violation, dark matter, and so on. In particular, it might be interesting 
to identify some of the KK right-handed neutrinos as the dark matter by introducing the 
bulk Dirac masses with an appropriate choice of the parameters. These subjects remain 
to be explored in future work. 
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A Lorentz spinors and gamma matrices 

In this work, the gamma matrices are taken as 

{T M ,T N } = 2 V MN = 2diag(+l, -1,-1, -1,-1), (A.l) 

rv = 7,= (-^ <r 4 = 7b=( 1 (A.2) 

where = (1,0",) and <r M = (1, — A 4-component spinor is written in terms of 
2-component spinors as 

* = (It) ■ (A.3) 



The Dirac and charge conjugates for \l/ are given by 

* = or a), * c = ^ T = fl^Q J > ( A - 4 ) 

where C5 is the charge conjugation matrix in five dimensions: C5 = 27 2 7°75. The anti- 
symmetric tensors are 

= = = (A-5) 



B S4 group 

The 54 group consists of all permutations of four objects and the order is 4! = 24. The 
5*4 is symmetry of the cube. The irreducible representations are two singlets 1 and 1', one 
doublet 2, and two triplets 3 and 3'. In this work, we adopt the following basis [18] 



«=l! J). p=( u Q $ cb.i; 



for the doublet 2 and 








V 











: 







!• H 


1 









1 


0) 




Vo 


1 





Q = 0-1, P = 1 (B.2) 
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for the triplet 3. All the group elements are given by the products of these two generators; 



(B-3) 



The tensor products 3x3 and 2x3 play an essential role to construct the Yukawa 
coupling matrices in symmetric phase. They are decomposed as 3x3 = 1 + 2 + 3 + 3' 
and 2x3 = 3 + 3'. Suppose that tpi, 0j (i = 1, 2, 3) and Xj (j — 1> 2) behave as 3 and 2 
under the basis of flB.2j) and fIB.lj) . respectively. Then it follows 



a i 


= Q\ 


a 2 


= Q 2 , 


a?, 


= PQ 2 P 2 , 




= Q 2 PQ 2 P 2 




= P, 


b 2 


= Q 2 P, 


b 3 


= QPQP 2 , 


h 


= Q 2 PQ 2 , 


Cl 


= P\ 


c 2 


= Q 2 P 2 , 


c 3 


= QPQ, 


c 4 


= Q 3 PQ, 


d\ 


= PQPQ 2 , 


d-2 


= PQP, 


d 3 


= Q 3 , 


d,\ 


= Q, 


ei 


= Q 2 PQ, 


e 2 


= PQ, 


e 3 


= Q 3 P 2 , 




= QP 2 , 


h 


= QPQ 2 , 


h 


= PQP 2 , 


h 


= P 2 Q, 


h 


= QP. 



%j) X 




^101 + W^202 + U 2 ljj 3 4> 3 
V>101 + UJ 2 1p2(p2 + ^303 



2 



^203 + ^302 \ / ^203 ~ ^302 
+ ( ^301 + ^103 + ^301-^103 | (B.l) 
^102 + ^201 / V ^102 - ^201 



3 3' 
(Xl + X2)tpl \ / (Xl ~ X2)^l 

X x ip = \ (u 2 xi + 0^X2)^2 + (w 2 Xi - WX2)^2 J • (B.5) 

(UXI + W 2 X2)^3 / V (WXI - ^ 2 X2)^3 



3 3' 
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